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ABSTRACT 



Context. We present a code, VeDyn, for modelling envelopes and atmospheres of A to F stars focusing on accurate treatment of 
convective processes. 

Aims. VeDyn implements the highly sophisticated convection model of Canuto and Dubovikov (1998) but is fast and handy enough 
to be used in practical astrophysical applications. 

Methods. We developed the HME envelope solver for this convection model furtheron to consistently model the envelope together 
with the stellar atmosphere. The synthesis code SynthV was extended to account for the resulting velocity structure. Finally, we tested 
our approach on atomic line bisectors. 

Results. It is shown that our synthetic line bisectors of HD 108642 bend towards the blue and are of a magnitude comparable to the 
observed ones. 

Conclusions. Even though this approach of modelling convection requires the solution of a coupled system of nonlinear differential 
equations, it is fast enough to be applicable to many of the investigation techniques relying on model atmospheres. 
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1. Introduction 

A physical mechanism most challenging to understand in F to 
early A stars is convection. It is one of the least understood phe- 
nomena taking place in these stars. Recently, several success- 
ful attempts to describe and to model convection in the enve- 
lope and the photosphere in more detail have been conducted. 
To be mentioned in this regard are for instance, numerical sim- 
ulations by Nordlu nd & S tein J2000b , Stein & N ordlun d d2003l 
Freytag & Steffen d2004b and Kochukhov et al. d20Q7T ), but also 
tests and applications of nonlocal convection models presented 
in Kupka ( TT^99l HME ) , Kupka & Muthsam d2007b . and Kupka 
& Montgomery J2002b . 

However, in practical applications like stellar spectrum anal- 
ysis or stellar evolution modelling the mixing length theory 
(Biermann 1 1948b is still most commonly used to model con- 
vection. Kupka d 19961 ) and Heiter et al. ([2002 ) implemented a 
local full spectrum turbulence model of Canuto and Mazzitelli 
dl991|1992| ) and Canuto, Goldman and Mazzitelli ( fl996l) in th e 
widely used model atmosphere code Atlas9 (Kurucz I1993I ). 
The complete Reynolds stress model (RSM) of Canuto and 
Dubovikov ( 1 19981 CD98) has been applied to envelopes of 
A-stars and also white dwarfs a few years ago (Kupka & 
Montgomery |2002l Montgomery & Kupka 2004). Their calcula- 
tions are in qualitative agreement with 2D numerical simulations 
of Freytag ( 119951 ) and the resulting vertical velocites are consis- 
tent with observerved values for micro- and macro turbulence. 

Extending this work, we investigated the possibility to ap- 
ply the CD98 convection model also to stellar atmospheres. The 
various analysis methods of stellar spectra connected to or rely- 
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ing on stellar atmosphere modelling mostly use the micro- and 
macroturbulence parameters to account for those contributions 
to the lines shapes believed to originate from turbulent veloci- 
ties. These however, permit only to model distortions symmetric 
with respect to the line center and are not inherent to the under- 
lying stellar model. 

In this publication we show our extension of the HME enve- 
lope solver to model stellar envelopes together with the stellar at- 
mosphere and a method to account for the derived velocity fields 
in spectrum synthesis. Following this approach, the structure and 
magnitude of the up- and downfiows is resulting from the model 
calculation and is then put into the derivation of spectral line 
profiles. Hence, the need for the parametrization of a micro- and 
marcroturbulence is expected to be considerably reduced. 

We also present an application of our models to the star 
HD 108642 which has already been investigated by Landstreet 
( 1 19981 ) and others concerning microturbulence, abundances and 
also bisectors. 

2. The models 

The convection model which we used is an extension of the 
CD98 approach, as described in Kupka and Montgomery ( 2002 ). 
The equations for the turbulent quantities are solved on a 
weighted mass grid. The weighting allows a non-equidistant 
spacing of the grid such that steep temperature gradients are re- 
solved. All other quantities are derived on the same weighted 
mass grid in the envelope and on a log(r ross ) grid in the atmo- 
sphere. The transition region is determined automatically ac- 
cording to p, P and T at each iteration step. This way we could 
consistently model the stellar atmosphere together with the en- 
velope with feasible computational efforts. For the models pre- 
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sented in this paper we assumed spherical symmetry. We aim to 
use at least two sources for the opacities which are OPAL opac- 
ities (1996 ) in the envelope, and opacitiy distribution functions 
in the atmosphere. At present, however, we restrict ourselves to 
the OPAL routines, since they are faster and basic code develop- 
ment is still ongoing. Depending on the maximum size of time 
steps and the necessary relaxation time one model requires 3 to 
30 hours to finish on a 2.3 GHz single CPU machine. 



3. Synthesis and bisectors 

To derive synthetic spectra from our models and velocity pro- 
files, we adopted the SynthV code of Vadim Tsymbal. Following 
the CD98 equations (38d) and (50d) to (50f) 
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we can derive typical velocities a>" p , <w^ own and the filling fac- 
tor cr (fractional area occupied by the updrafts) from the mean 

vertical turbulent velocity to 1 and the skewness S w . 




Fig. 1. Velocities (-w" p - dashed, of 
derived from our model. 
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Figure Q] compares the absolut values of the up- and 
downdraft velocities we calculated. Note that at tau = 2/3 
(log(r loss ) = -0.176) these are in the range of typical values for 
the microturbulent velocity one would expect for an A type star. 
Contrary to CD98 our upstream velocities are negative, since in 
our frame of reference up means moving towards the observer. 

The synthesis code SynthV of V. Tsymbal has been de- 
veloped further on to allow to account for the upward and 
downward flow by means of separate radial velocity profiles. 
Their contributions to the total line profile are then summed up 
weighted according to the filling factor (Figure |2). 
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Fig. 2. Filling factor in our model of HD 108642. 



4. HD 108642 

The Am star HD 10842 has been observed several times over the 
past years concerning its peculiar nature, abundance pattern and 
also (Landstreet l 19981 1 atmospheric velocity fields. The observa- 
tions we used for this investigation were done by Landstreet and 
Kupka in April 2001 at the CFHT They observed the star at a 
very high resolution of R = 120000 in several chunks of approx- 
imately 100 A in range. We found the ranges 5491 - 5572 A and 
6102 - 6191 A most suitable. These contain enough unblended 
lines at a reasonable noise level S/N > 140 and 150. 

As starting point for our model for the envelope and atmo- 
sphere of HD 108642 we took evolutionary envelopes as well 
as polytropes. We found our code to be insensitive to the dif- 
ferences of these types of envelope models. Extension into the 
atmosphere was performed with the gray approximation. For ra- 
diation transport we used the OPAL opacities and equation of 
state Rogers et al. 1996) as input data. The only tuneable pa- 
rameters we had in our model of HD 108642 were its effective 
temperature of 8100K, the stars surface gravity of 4.1 and we 
accounted for its higher metallicity by using opal tables for Z 
= 0.06. All other parameters and numerical constants are as de- 
scribed in Kupka and Montgomery (2002). To do the line syn- 
thesis we adopted the SynthV code of Tsymbal. Its approach of 
solving the equation of radiative transfer for each atmospheric 
layer separately naturally allows to account for depth depen- 
dent velocity fields. We performed the radiative transfer twice 
for each layer, for the upwards and for the downwards moving 
matter. These two contributions were then summed up weighted 
according to the filling factor and the absorption of the individ- 
ual layer added to the total line absorbion. 

FigureQ] shows the mean velocities of upwards and down- 
wards moving material resulting from our calculations. The val- 
ues of a>" p and w down taken at T ross = 2/3 compare reasonably 
to the measurement of v m ; c = 4.0kms~ 1 by Landstreet (I19981 ). A 
comparison of the observed bisector of the CrII line at 4616 A 
(Landstreet 119981 ) and the synthetic bisector derived from the 
VeDyn model can be seen in Figure|3] Clearly the magnitude 
of the calculated bisector is of the same order as the observed 
one. Also the tilt to the blue is reproduced. However, its shape 
is quite deviating. This, we believe, is at least in part due to the 
fact that we could not include the upper part of the atmosphere 
(log(T, OSJ ) < -4.0), because of numerical problems. Furthermore 
figure [3] shows the bisectors of two weak lines, Col 5342.7 and 
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Fel 5241.9; an encouraging result of our modelling procedure. 
However, we currently do not have any observations to compare 
those synthetic bisectors to. 




-0.04 -0.03 -0.02 -0.01 0.01 



wave - wave_center [A] 

Fig. 3. The solid line shows the observed bisector of CrII as 
measured by Landstreet 1998, the dashed line denotes the syn- 
thesized bisector. The dotdashed and dotted line represent the 
bisectors of Col 5342.7 and Fel 5241.9 respectively; both weak 
lines. 
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5. Conclusions 

We consistently modeled turbulent convection in the lower at- 
mosphere and envelope of an Am star. Following the model of 
Canuto and Dubovikov (1998 ) and the approach of Kupka & 
Montgomery (2002) we solved the equations for the turbulent 
quantities on a weighted mass grid. Other quantities are com- 
puted on the same weighted mass grid in the envelope and on a 
log(r ross ) grid in the stellar atmosphere. Line synthesis has been 
performed with a modified version of the SynthV code of Vadim 
Tsymbal, where we account for the depth dependent mean up- 
and downstream velocities and the filling factor. 

The resulting line shapes of HD 108642 have been compared 
to the observed ones. The same has been done for the typical 
turbulent velocities we derive at optical depths in the vicinity of 
r ross = 2/3. They are in the range of the microturbulent velocity 
measured by Landstreet ( 1998]). 

We also calculated bisectors from our synthetic lineprofiles. 
They are tilted to the blue by a magnitude comparable to the ob- 
served ones. However, the shape of the synthetic bisector differs 
from the observation. We think this results from the upper part 
of the atmosphere missing in our model. Currently more objects 
are investigated to cover the region of A to F type main sequence 
stars. 

Our modelling and synthesis approach is applicable to inves- 
tigation methods relying on model atmospheres. Still, the fre- 
quency dependent radiation transport has to be implemented in 
our models. 
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